Theory of concentration dependence in drag reduction by polymers and of the MDR 
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A simple model of the effect of polymer concentration on the amount of drag reduction in tur- 
bulence is presented, simulated and analyzed. The qualitative phase diagram of drag coefficient vs. 
Reynolds number (Re) is recaptured in this model, including the theoretically elusive onset of drag 
reduction and the Maximum Drag Reduction (MDR) asymptote. The Re-dependent drag and the 
MDR are analytically explained, and the dependence of the amount of drag on material parameters 
is rationalized. 



"Drag reduction" refers to the intriguing phenomenon 
when the addition of few tens of parts per million (by 
weight) of long-chain polymers to turbulent fluids can 
bring about a reduction of the friction drag by up to 
80% 0, 0, Q The phenomenon is well documented 
since Toms discovered it accidentally in 1946 while study- 
ing the degradation of polymers. The pioneering work 
of Virk 0, had systematized and organized a huge 
amount of experimental information, but the fundamten- 
tal mechanism for the phenomenon has remained under 



MDR 



debate for a long time 
many of the numerical 




. All the experimental and 
investigations of drag 



reduction focused on channel and pipe geometries; re- 
cently however it had been discovered by numerical sim- 
ulations ^3 of model equations of viscoelastic flows (like 
the FENE-P model) that drag reduction appears also in 
homogeneous and isotropic turbulence when seeded with 
polymers. This brought about a new focus to the search 
for the mechanism for drag reduction, since the analy- 
sis of model equations without wall effects should suffice 
to uncover a mechanism. Indeed, in a recent paper [llj 
the FENE-P equations were simplified further to a shell 
model of viscoelastic flow which was shown to exhibit 
drag reduction whose mechanism could be fully explored 
analytically. In this Letter we present additional crucial 
progress where we demonstrate and explain two of the 
most prominent (and least understood) characteristics of 
drag reduction, i.e. the onset [as a function of Reynolds 
number (Re)] and the Maximal Drag Reduction (MDR) 
asymptote. 

To set up the issues we reproduce in Fig.l a typical 
experimental figure from Ref. @ which refers to the de- 
pendence of the friction (or drag) coefficient in pipe flows 
on Re. For a pipe of radius R and length L, with Ap, p 
and U being the pressure drop across L, the fluid density 
and the mean velocity over a section, the drag coefficient 
/ reads 




Prandtl - Karman 
Law 



Ap R 



(Pipe flow) 



(1) 
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FIG. 1: Drag reduction in Prandtl-Karman coordinates Q. 
As a function of Re the drag exhibits a (concentration inde- 
pendent) transition to drag reduction. The amount of drag 
reduction depends on the concentration until the asymptote 
denoted by MDR is reached. The Prandtl-Karman law is the 
Re-dependent drag of the neat fluid. The numbers indicate 
concentrations of the polymer additive in wppm. 



"Drag reduction" is tantamount to, say, an increase in the 
throughput U for a fixed pressure drop Ap when polymer 
is added to the working fluid. In Fig. 1 one sees that for 
low Re there is no drag reduction: the drag coefficient 
of pure water is unchanged by the addition of small con- 
centration of polymers. Then there is a sharp onset of 
drag reduction at a value of Re that does not depend on 
the concentration. From this point on the amount of drag 
reduction depends both on the concentration of the poly- 
mer and on Re. It was shown by Virk however that the 
amount of drag reduction asymptotes to an apparently 
universal curve that cannot be exceeded by increasing 
the concentration further. This asymptote is referred to 
as the MDR, and was claimed to be insensitive to the 
nature of the polymer used in the experiments. In spite 
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of the ample experimental evidence, both the onset and 
the existence of the MDR have not been theoretically 
understood. In this Letter we wish to close this gap. 

Our strategy is to explore simulationally and analyti- 
cally simplified models of viscoelastic flows which in spite 
of the simplification still represent the robust properties 
that we are after. As is well known, viscoelastic flows 
are represented well by hydrodynamic equations in which 
the effect of the polymer enters in the form of a "confor- 
mation tensor" R(r, t) which stems from the ensemble 
average of the diadic product of the end-to-end distance 
of the polymer chains. Flexibility and finite extendabil- 
ity of the polymer chains are reflected by the relaxation 
time r and the Peterlin function P(r, t) which appear in 
the equation of motion for R. 
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In these equations p 2 m and p\ refer to the maximal and 
the equilibrium values of the trace i? 77 . Since in most 
applications p rn 3> po the Peterlin function can be also 
written approximately as P(r,t) w 1/(1 — ai? 77 ) where 
a = p~ 2 . In its turn the conformation tensor appears in 
the equations for fluid velocity u(r,t) as an additional 
stress tensor: 



T(r,t) = 



P(r,t) 



v s V 2 u 



v-t + f 



2 fl(r,t)-l 

Po 



■(4) 
(5) 



Here v s is the viscosity of the neat fluid, F is the forcing 
and v p is a viscosity parameter which is related to the 
concentration of the polymer, i.e. v v jv s ~ c where c is 
the volume fraction of the polymer. Note that the ten- 
sor field can be rescaled to get rid of the parameter a in 
the Peterlin function, R a p — aR a p with the only conse- 
quence of rescaling the parameter po accordingly. These 
equations were simulated on the computer in a channel or 
pipe geometry, reproducing faithfully the characteristics 
of drag reduction in experiments 0, 0, 01 ■ It should be 
pointed out however that even for present day computers 
simulating these equations is quite tasking. We therefore 
simplify the model further. 

In developing a simple model we are led by the fol- 
lowing ideas. First, it should be pointed out that all 
the nonlinear terms involving the tensor field i?(r, t) can 
be reproduced by writing an equation of motion for a 
vector field B(r,t), and interpreting R a p as the diadic 
product B a Bp. The relaxation terms with the Peterlin 
function are not automatically reproduced this way, and 
we need to add them by hand. Second, we should keep 
in mind that the above equations exhibit a generalized 



energy which is the sum of the fluid kinetic energy and 
the polymer free energy. Led by these consideration we 
write the following shell model [ll| , which we refer to as 
the SabraP model: 

^ = U n (u,u)- l -^P(B)$ n (B,B)-v s k 2 n u n + F n , 
at o 6 t 

dB iil 

= -$ n (u,B)--$ n (B,u)--P{B)B n -v B k 2 l B n , 
at o o t 



P(B) = 



1 



(6) 



In these equations u n and B n stand for the Fourier am- 
plitudes u(k n ) and B(k n ) of the two respective vector 
fields, but as usual in shell model we take n = 0, 1, 2, . . . 
and the wavevectors are limited to the set k n = 2". The 
nonlinear interaction terms take the explicit form 



$„(u, B) = k n [(1 - b)u n+2 B* +1 + (2 + b)u* n+1 B n+2 
+fc n _i (26 + lX^Sn+i - (1 - b)u n+ iB*_ 1 
+k n -2 (2 + 6)u„_iS„_2 + (26 + l)w n _ 2 -Bri-l 



(7) 



with an obvious simplification for <&«(«, u) and <!>„(£>, B). 
Here b is a parameter taken below to be —0.2. In accor- 
dance with the generalized energy of the FENE-P model, 
also our shell model has the total energy 
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The second term in the generalized energy contributes 
to the dissipation a positive definite term of the form 
(v p /t 2 )P 2 (B) J2 n \B„\ 2 . With v p = the first of Eqs. 
reduces to the well-studied Sabra model of Newtonian 
turbulence 01- ^ s m the FENE-P case we consider 
v v jv s to be c. All the simulations below are performed 
with a constant rate of energy input, choosing F n = <j>/u n 
for n = 0, 1 and zero otherwise. 

In ^l| it was shown that this shell model exhibits drag 
reduction, and the mechanism for the phenomenon was 
elucidated. The basic phenomenon is exhibited well by 
the spectra of the u n and B n fields which are presented 
at one value of the parameters in Fig. The spectra 
for the Sabra model (dashed line) and the SabraP model 
(line) are compared for the same amount of power input 
per unit time. The discussion |ll| of the spectra revolves 
around the typical Lumley scale k c which is determined 
by the condition 



u{k^)k c 



(9) 



For k n 3> k c the decay time r becomes irrelevant for the 
dynamics of B n . The nonlinear interaction between u n 
and B n at these scales results in both of them having the 
same spectral exponent which is also the same as that 
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c=0.5 - 
c=0. 



FIG. 2: Power spectra of the SabraP model (line) and the 
Sabra model (dashed line) for <f> = 0.001, v s = 10~ 6 and 
r = 0.4. The dashed line with symbols represents the power 
spectrum of the B„ field. 



of the Sabra model. The amplitude of the u n spectrum 
is however smaller in the SabraP model compared to the 
Sabra case, since the B n field adds to the dissipation. 
On the other hand, for k n <C k c the B n field is exponen- 
tially suppressed by its decay due to r, and the spectral 
exponent of u n is again as in the Sabra model. Drag 
reduction comes about due to the interactions at length 
scales of the order of k c which force a strong tilt in the u n 
spectrum there, causing it to cross the Sabra spectrum, 
leading to an increase in the amplitude of the energy con- 
taining scale. This is why the kinetic energy is increasing 
for the same amount of power input, and hence drag re- 
duction. Note that a very similar spectral cross-over had 
been documented also for the FENE-P model in channel 
flow simulations @ 

The qualitative phenomena that we are about to ex- 
plain in this Letter are demonstrated in the simulational 
results presented in Fig. Here we show the drag co- 
efficient / as a function of Re for the Sabra and for the 
SabraP models for various values of the concentration. 
The drag coefficient is computed in analogy to Eq. 
as 



f = 



(e, 



3/2 

121 fco 



(Our model) (10) 



We observe all the phenomena discovered by Virk: (i) For 
the model of the neat fluid the drag has a laminar branch 
and a turbulent branch, with a sharp transition between 
them, (ii) For the model of the viscoelastic flow in the 
laminar region there is no drag reduction; the laminar 
branch is not changed by the addition of polymer with 
any concentration, (iii) Drag reduction has an onset that 
is independent of the concentration of the polymer, (iv) 
As the concentration increases the amount of drag reduc- 
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FIG. 3: Upper panel: Drag as a function of log 10 (Re) in- 
cluding the laminar and the turbulent regimes. Both regimes 
agree with the Eas. (|1 (|12|l . Lower panel: Blow up of the 
turbulent regime. In both panels the upper straight line indi- 
cates the drag of the neat fluid, whereas the MDR is seen as 
the convergence of the drag data for large concentrations. 



tion increases, but (v) there exists an asymptote which 
is not exceeded when the concentration is increased. In 
other words, our simple model appears to reproduce ex- 
tremely well the phenomena that were uncovered in so 
many experiments as summarized by Virk. 

Next we explain all these observations. First we ra- 
tionalize the Re-dependence of the friction factor in the 
Sabra model of the neat fluid. For low Re the nonlinear 
terms $> n (u, u) are negligible compared to the viscous 
term. Forcing only on the largest scale kg we can evalu- 
ate, 



Re" 1 Re small . (11) 



For large Re we have the exact result |l2| that the third 

to} 

order correlation function Sn = Q(u n -iu n u^ +1 ) = 
Ce/k n with e being the mean energy flux and C a known 
constant (the analog of the 4/5th law for Navier-Stokes 
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turbulence). We therefore expect the friction factor 
to tend to a constant value for large Re (up to terms 
~ logRe), 

/ - Re° Re large . (12) 

Eqs. jnjl an d (which are the analogs of the Prandtl- 
Karman law for pipe flows) are well borne out by the 
data in Fig. [21 for the model of the neat fluid. The lam- 
inar branch, which is exponential in these coordinates 
(/ ~ exp{— log(Re)}), is unaffected by changing the con- 
centration c. The transition between the two branches is 
expected when turbulence sets in, i.e. for Re such that 
the dissipative terms just begin to be overwhelmed by 
the nonlinear interactions. Thus point (i) is understood. 
Note that similar arguments will hold for the FENE-P 
equations in homogeneous flows. Points (ii) and (iii) are 
explained as follows; we said above that drag reduction 
comes about due to the interaction between the two dy- 
namical fields at scale of the order of k c . Clearly, as long 
as k c exceeds the dissipative scale kd of the the u n field, 
no interaction between the two field can be of any signif- 
icance. Since kd is of the order of kd ~ kg -Re 3 / 4 , we can 
expect a concentration independent onset of drag reduc- 
tion when k c w kd- Using u n ~ uo(^n/fco) -1 ^ 3 , k c can be 
estimated as k c ~ ko(TUoko)~ 3 ^ 2 , and we end up with a 
prediction for the onset of drag reduction when 

Re w (ruofco) -2 Onset of drag reduction . (13) 

This prediction is well borne out by our simulations (due 
to the space constraint we do not display simulations 
at different values of r and fco). Again we point out 
that similar arguments can be presented for the FENE-P 
model as well. 

Point (iv) is obvious - when the concentration in- 
creases the mechanism discovered in |ll| comes into play 
What remains to explain is the asymptotic MDR. This 
also follows directly from the analysis of the equations. 
Consider Eqs. © for two values of the parameter v p , 
<C fp 2 \ with y 2 = Up jv^ . Rescaling B n according 
to B n = yB n , we see that the Peterlin function tends to 
unity when y — > 0, 

P ^) = - * ,5 tz ->1 , wheny^O. (14) 

When P{B) w 1 the dynamical equation for B n is in- 
dependent of of y due to its linearity, whereas the u n 
equation remains independent due to the rescaling: 

■ (2) ■ (1) 

l-E-P(B)$ n (B,B) - i-H-P(B)$ n (S n> B„) , (15) 
3 t 3 r 



Thus increasing the concentration brings to the dynam- 
ical equations to an asymptotic concentration invariant 
form and therefore to an asymptotic MDR. For the last 
time we remark that similar rescalings are also available 
in the FENE-P equations, making the points discussed 
here quite general for any sensible model of viscoelastic 
flow. 

In summary, we have presented a simple model of drag 
reduction for which the observed characteristics can be 
explained on the basis of the equations of motion. It 
remains to go back to channel and pipe simulations of the 
FENE-P equations to demonstrate that the discussion 
presented above includes the main phenomena observed 
also there. 
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